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In this paper we study two-dimensional inside-out Eaton lens theoretically and numerically. With 
the help of the WKB approximation, we investigate the finite- wavelength effect analytically and 
demonstrate one necessary condition for perfect imaging by the lens. Furthermore, we present one 
example design and test its performance by using full wave Maxwell solvers. 



More than half a century ago, it was proposed that gradient index lenses [1], such as the Maxwell fish-eye lens 
[2], the Luneburg lens [3] and the Eaton lens can be free of geometrical aberrations and form perfect images, at 
least at the geometrical-optics level (see Reference @ and [g^ for more details). Taking the inside-out Eaton lens as 
an example, its refractive index n(r) equals ^J{2 — r)/r for 1 < r < 2 and 1 otherwise, with r being the radius (see 
Figure 1). We can analytically prove that light rays emitted from a source at position tq, with ro < 1, will be focused 
exactly at position — tq [6.]. Therefore, like Maxwell's fish-eye lens, perfect imaging can be obtained by an inside-out 
Eaton lens, and both the source and the image are inside the optical instruments. 

Recently, a renaissance of scientific interests appears in these gradient index lenses 04121 , partially because of 
the developments of metamaterials [13., 1^ and transformation optics [l5|-[l8[. Metamaterials are manmade media 
whose effective permittivities and permeabilities are determined by their deeply subwavelength structures as well 
as their constituent materials [l^. For instance, integrating split ring resonators with metallic rods will result in 
a metamaterial with an effectively negative index of refraction [19]. Transformation optics, on the other hand, is 
inspired by an intriguing property of Maxwell's equations, i.e., their form is invariant under arbitrary coordinate 
transformations, assuming the field quantities and the material properties are transformed accordingly jlBI, fl6|. 

In this paper, we design a two-dimensional dielectric inside-out Eaton lens consisting of metallic wires in a homoge- 
neous background medium with positive permittivity. To investigate the effect of the finite wavelength, we employ the 
WKB approach to solve the wave equations [111, and obtain one necessary condition for perfect imaging. Finally 
we test our design with full wave Maxwell solvers. 

The original Eaton lens has a spherical geometry and an index of refraction as 



n(r) = ^'-l, (1) 

with r being the radius [i*]. Since it possesses spherically symmetry, the entire ray trajectory lies in a plane which is 
orthogonal to a conserved angular momentum L [6] (See Appendix A). Consequently, we can simplify it to be a two- 
dimensional cylinder with similar refractive index. We further assume that the cylinder lies in the xy plane, as well as 
the propagation plane of the light ray. Moreover, we only consider TM-polarized light where the magnetic-field vector 
H points in the z direction. By sacrificing the impedance matching, we can further assume that the two-dimensional 
cylinder is purely electrical, i.e. = with its permittivity e(r) given by n^{r). 

A detailed ray-optics description of the inside-out Eaton lens can be found in Ref. [6] , and it is repeated in Appendix 
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FIG. 1: Inside-out Eaton Lens (1 < r < 2). The source is located at ro = 0.5. Light rays (blue) are described by Hamilton's 
equation (A. 2). 
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FIG. 2: (a) Dependence of the function s(r) on the mode order n. Here the wavelength is 0.3. (b) Effect of the wavelength A 
on the phase factor 7n. 



A for the convenience of readers. We therefore only consider the wave interpretation here. Starting from the following 
wave equation 



e[r) e[r) 

with = dx^x -\- dyBy and H as H^^ this equation can be reformulated in cylindrical coordinates as 



r e dr J dr dO'^ ^ ' 



(2) 



(3) 



with ho = uo/c being the wave number in free space. We now assume that the magnetic field i^(r, 0) can be expanded 
as Xln /n(^)e*^^\/e/r, where the functions fn satisfy 



n 



In the region r < 1 where e(r) = 1, we have 
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and express the solutions generally as 



fn{r) = \anHi^\kor) + bnH^^\kor) 
rder ] 

1 < r < 2, we can rewrite Equation (|4]) as 
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where Hn^ and Hn^ are the n-th order Hankel functions of the first and second kind, respectively. In the region 



with 



s(r) = 



fn + s{r)klU = 0, 
2-r 



r + 1 



(7) 
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A few examples are plotted in Figure (2). For a modest n, s{r) generally monotonically decreases from a positive 
value to negative infinity. It is however always negative when n is large enough. 

We employ the WKB approximation, developed by Wentzel, Kramers and Brillouin in 1926, to analytically solve 
Equation ([7]) with a modest n [20|, |2l|. More specifically, we assume that fn{r) has the form A^e^^^^^'^^ for positive 
5(r), and r(r) can be further expanded in terms of /cq. 



r{r) = ro(r) + i^n(r) + T2^2(0 

/Co /VQ 



(9) 
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FIG. 3: (a) A schematic of the design, (b) The radius of the metalhc wire as a function of its position. Here the permittivity 
of the metal is Cm — —0.6. 



Similar arguments also hold by t' {r) as well as t" {r). By collecting the leading-order terms, it is found that 



dr- 



dr 

Consequently the first order solution can be expressed as 



dr 2r/, 
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with Vn standing for the turning point where s{r) = 0. The first term on the right hand side corresponds to an 
out-going wave because its phase increases with distance, while the second term corresponds to an in-coming wave. 
Similar procedures can be applied to a negative s(r) by assuming /n(r) = Bne~^^^^^\ and the resultant first-order 
approximation is given by 



/n(r) 



|s(r)|V4 



exp 



(12) 



Here only the solution that is exponentially decaying in the r direction is included. Furthermore, s{r) ~ (r^ — r) 
approaches zero linearly in the vicinity of the turning point Vn- The solution therefore can be approximated as 
Ai(— /cq^^s), with Ai being the Airy function poj. When /cq is large enough, we can asymptotically match Equation 
(pT]) and ([12]) around the turning point, and finally achieve = iA~. As a direct result, in the region where s{r) is 
positive, we have 



fn{r) 



5(r)V4 



exp 
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which implies that an out-going wave will be totally reflected around the turning point r^, accompanied with a phase 
variation of 7r/2. 

We now assume Equation (p!3|) can be extended to the region where r is slightly smaller than 1. Additionally, its 
phase factor can be approximated as 



^0 j \/s{r')dr' = ko j ^/%dr + /cq / \/ s{r') - ^/s^ dr ^ kor + 7n, 
where 7^ does not depend on r, and 



^0"^0 — ^0 



n2 - 1/4 



being the coefficient shown in Equation In the vicinity of r = 1, fn{r) can be then rewrite as 

gi(fcor+7„) _|_ g-i(fcor+7„+7r/2) 



fn{r) ~ 



(14) 
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FIG. 4: Full-wave simulations of the design, with the excitation source placed at different locations, ro equals (a) 0.2, (b) 0.3, 
(c) 0.4, and (d) 0.5. The wavelength is fixed at 0.3, and the permittivity of metal is Cm = —0.6 + O.Oli. The amplitude of the 
magnetic field |i^| is plotted. 



Coincidentally, a source at position ro with ro < 1^ i.e. inside the Eaton lens, generates a radiation such as 



ff(^'(fco|r - rol) = -T E Jn{kor<)Hi'\kor>)e'"t 



(17) 



where r< is the smaller of r and ro , and r> is the larger of r and ro [201, l22| . The total magnetic field between ro and 
1 is therefore given by 



(18) 



with Cn representing the amplitude of the reflected n-th order wave. When ko is large enough, the large argument 
approximations of the Hankel functions and the Bessel function lead to 



koirJrr^ 



cos(/coro - Pn 



(19) 



with /3n = (2n + l)7r/4. Again, we asymptotically match the above equation with /n(^)/v^ of Equation (p!6|) . and 
finally achieve 



Cn 



i[(n+l)7r+27„ 



(-1) 



n+1 -2i7„ 



7 ^r, 



25, 



1/4 



cos(/coro - /3n). 



(20) 



As mentioned at the beginning, an inside-out Eaton lens will convert light rays emitted from a source at position 
To, with ro < 1, to the opposite location — tq and result in a perfect image. To extend this property to waves, it is 
required that 



Cn = (-1) 



(21) 



with (j) being the phase difference between the source and image. This relation can be easily obtained by time reversing 
the source radiation process [22]. Comparing Equation (|2T]) with Equation (|2Q]) . we obtain the following necessary 
condition for perfect imaging by an inside-out Eaton lens: e*^^'^ must be constant and independent of the mode order 
n. The function e*^^'^ can then be employed to partially evaluate the performance of a lens. For instance, we calculate 
the phase factor 7^ by using Equation (fT4|) . and the results are shown in Figure (2b). Evidently, 7^ are almost 
constant when A is close to zero (the ray-optics region), while they vary strongly for a considerable wavelength. 

In the following, we present a conceptual design consisting of a vacuum (background) and ideal metalhc wires 
(a realistic device is currently under construction). First, the Maxwell- Garnett formula is used to approximate the 
effective permittivity of the composite (See Appendix B) 



(l-/) + (l + /)e^ 
(l-/)e^ + (l + /)' 



(22) 
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with / being the fihing fraction of the wires and Cm being the permittivity of the ideal metal [22|, |23| . Since eg should 
be equal to the permittivity of the Eaton lens, (2 — r)/r, the filling fraction is then given by 

/(r) = (r-l)l±^. (23) 

Notice that the first-order surface mode will be excited when e^n = — 1 [23]. The designed lens, shown in Figure (3a), 
consists of 8 layers of metallic wires with different radii. The above equation is then used to determine the positions 
as well as the radii of these wires. The thickness of the p-th layer a^+i — ap is set to be ap7r/30, with ap and a^+i 
being the inner and outer radius of the layer, respectively. Consequently we have 

ap = ai{1^7r/30y-\ (24) 

where ai is assumed to be 0.95. Furthermore the following relation 

is employed to calculate the radius of the wire in the p-th layer. In the current design, = —0.6 + iS with 6 being 
very small, which leads to the wire radii shown in Figure (3b). Evidently, the wire radius increases rapidly with the 
increasing of r, with a minimum of 0.007 and a maximum of 0.06. 

We employ a finite-element full- wave Maxwell solver to verify the above design [23], and plot the results with 
6 = 0.01 in Figure (4). Notice that we always set A = 0.3, a value much smaller than the lens size while big enough 
to ensure the validity of Equation (|22]) . Clearly, an image is always observed at the location opposite to the source. 
To evaluate the quality of the image, we define 



Hi-ro) 



(26) 



the ratio of the intensity of the image to that of the source. The bigger the 77, the better the lens is. In Figure (4), 
the source location tq is gradually increased from 0.2 to 0.5, with an increment of 0.1. The corresponding 77 is found 
to be 0.21, 0.27, 0.15 and 0.15, respectively. Along the azimuthal direction, we observe a few intensity minima, and 
the total number of these minima depends on the source location tq. This phenomenon is very likely induced by 
the finite wavelength A. One direct consequence is that different modes radiated from the source have quite different 
amplitudes, as given by Jn{koro) of Equation (pT]) . For example, the zeroth-order and third-order modes dominate 
the radiation when ro = 0.2, while the ninth-order mode is the strongest one when ro = 0.5. We further investigate 
the infiuence of the metallic absorption in Figure (5), by setting ro = 0.2 and A = 0.3. Two different 0.1 and 0.5, 
are considered, and the corresponding intensity ratio r] is found to be 0.11 and 0.06, respectively. Evidently, although 
the image quality is degraded with the increasing of the metallic absorption, the basic function of the Eaton lens, 
i.e. forming an image, still survives. We can partially interpret it with the fact that all the higher-order modes, such 
as n > 25, are totally refiected around r = 1. Only the low-order modes can propagate into the lens and hence be 
absorbed by the metal. 

To summarize, we have studied two-dimensional inside-out Eaton lenses. The corresponding full wave equation is 
analytically solved with the help of the WKB approximation. One necessary condition for perfect imaging is further 
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found, ie, e*^'^'^ must be independent of the mode order n. Furthermore, a general design procedure of the lens, 
based on effective medium theory, is developed. We present one example consisting of metal wires with different radii, 
and further verify the design with a full wave Maxwell solver. Its dependence on source location as well as metallic 
absorption is also investigated. 

This work was supported in part by the Penn State MRSEC under NSF grant no. DMR 0213623. 



Appendix A: A ray-optics theory of the Eaton lens 

A detailed description regarding the ray-optics theory of the Eaton lens can be found in Reference [6^. We briefly 
repeat it here for the reader's convenience. 

In geometric optics, there are two different but equivalent ways to describe the trajectory of a light ray. The first 
one is the Newtonian Euler-Lagrange equation 



d'^T _ Vn2(r) 
d^ 



where n is the refractive index and the parameter ^ is given by d^ = dr/n. We can interpret the above equation by 
using Newton's law, ma = — Vt/, for a mechanical particle with unit mass moving in "time" ^ under the influence of 
potential U = — n^/2 + with E being an arbitrary constant. The second way is based on Hamilton's equation 



dr ck dk ck^ , . , , 

dt nk dt 



with k being the wave vector and c being the speed of light in free space. Notice that by treating frequency uj = ck/n 
as the Hamiltonian, the above equation resembles the standard form of Hamilton's equation. 
We can define an angular momentum as 

^ dr n , / * X 

L = rx- = -rxk, (A.3) 

which leads to 

dL d^Y 1 „ 9 din? r x r , , 

-— = YX -— = -YX Vn^ — 0, (A.4) 

when the refractive-index profile n(r) is spherically symmetric. The above equation suggests that the angular momen- 
tum L is conserved. Hence, a family of light rays propagating in the xy plane at the beginning will always stay in the 
same plane. This fact implies that a two-dimensional Eaton lens with similar refractive-index profile n(r) functions 
identically to the three-dimensional one. 

To solve the two-dimensional Newtonian Euler-Lagrange equation, it is convenient to introduce the complex number 
z = X -\-iy^ and further reformulate the equation as 

d'^z z dr? 11 



by substituting the refractive index of the Eaton lens n(r) = \/ (2 — r)/r. The solution, following Equation (6.13) 
and (6.14) of Reference Q, can be expressed as 

z = e^" [cos(2^') + ^sin7sin(2^0 +C0S7] , dS, = 2\z\dS,\ (A.6) 

which describes displaced ellipses rotated by the angle a. 



Appendix B: Maxwell- Gar net t Formula 

Consider a two-component mixture composed of inclusions embedded in an otherwise homogeneous matrix, where 
Cm and Cd are their respective dielectric functions. The average electric field (E) over one unit area surrounding the 
point X is defined as 

(E(x)) = J /^E(x')dx' = /(E„(x)) + (1 - /){E,(x)), (B.l) 
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with / being the volume fraction of inclusion. A similar expression can be obtained for the average polarization 

{P(x)) = /(P™(x)) + (1 - /){P,(x)). (B.2) 

We further assume that the fohowing constitutive relations are vahd 

(P™(x)) = eo(em - l)(E™(x)), {Pd(x)) = eo(ed - l)(Ed(x)), (B.3) 

and the average permittivity tensor of the composite medium is defined by 

(P(x))=eo(ee-I)-(E(x)). (B.4) 

Combining the above equations we can obtain the effective permittivity eg. Clearly the resultant eg depends on the 
relationship between (E^(x)) and (E(i(x)) |25i] . 

We now assume that the inclusion has the shape of a cyhnder, and its radius is far smaller than the wavelength so 
that its optical properties can be well described by the electrostatic equation 

V • (e(r)^) = 0. (B.5) 

By matching the boundary conditions we can prove that ^m/^o = "^^d/i^d + e^), where (j)rn is the total potential 
inside the cylinder when the external electric field — V^o is homogeneous. This relation is further used to obtain the 
electric field [26] . It is finally found that the average permittivity is scalar and can be expressed as 

consistent with the Maxwell- Garnett dielectric function. Equivalently we can express the filling fraction as 

/(r) = (b.7) 
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